Global Phase Diagram of z/ = 2 Quantum Hall Bilayers in Tilted Magnetic Field 
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We consider a bilayer quantum Hall system at total filling fraction v = 2 in tilted magnetic field 
allowing for charge imbalance as well as tunneling between the two layers. Using an "unrestricted 
Hartree Fock," previously discussed by Burkov and MacDonald (Phys Rev B 66 115323 2002), we 
examine the zero temperature global phase diagrams that would be accessed experimentally by 
changing the in-plane field and the bias voltage between the layers while keeping the tunneling be- 
tween the two layers fixed. In accordance with previous work, we find symmetric and ferromagnetic 
phases as well as a first order transition between two canted phases with spontaneously broken (7(1) 
symmetry. We find that these two canted phases are topologically connected in the phase diagram 
and, reminiscent of a first order liquid-gas transition, the first order transition line between these 
two phases ends in a quantum critical point. We develop a physical picture of these two phases and 
describe in detail the physics of the transition. 



I. INTRODUCTION 

Over the last fifteen years, one of the most exciting 
frontiers in two dimensional electron physics has been the 
study of quantum Hall bilayers^. Stacking two quantum 
Hall systems a small distance away from each other serves 
two main purposes: First of all, it is a way of creating a 
multicomponent system in which the layer index of the 
electron plays the role of an isospin. Secondly, it is a 
step in adding another dimension to the traditionally two 
dimensional quantum Hall medium. 

One of the most striking phenomena that illustrates 
the richness of the physics of bilayer quantum Hall sys- 
tems is the commensurate-incommensurate transition ob- 
served in systems with filling fraction i/ = 1 subjected to 
tilted magnetic fieldsi^SiSi^i^. For typical parameters in 
these systems, the spin degrees of freedom are frozen out 
and the only important discrete degree of freedom is the 
layer index (isospin). Thus only two processes are im- 
portant: interlayer tunneling and Coulomb interactions. 
When such a bilayer system is subjected to tilted mag- 
netic fields, the electrons tunneling between the layers en- 
close flux quanta. In the presence of a sufficiently strong 
in-plane magnetic field it may be favorable for the system 
to forgo tunneling in order to escape destructive interfer- 
ence induced by the in-plane field. This creates an op- 
portunity for a phase transition to occur between a phase 
in which tunneling plays the central role and a phase in 
which tunneling is effectively zero. 

Recently, it has been showr^il that a somewhat differ- 
ent phase transition can occur in bilayer systems with a 
total filling fraction v = 2. The = 2 bilayer system is 
enriched not only by the presence of two electrons per 
flux quantum, but most importantly by the role of real 
spin. In contrast to = 1 bilayers where the real spin 
degrees of freedom are frozen out, in v ~ 2 bilayers, spin 
degrees of freedom are entangled with the isospin degrees 
of freedom as a result of Fermi statistics (Pauli exclusion 
principle, in other words). Thus, the v = 2 bilayers ex- 
hibit a richer phase diagram even in perpendicular field. 



and, as we will see later, an even more intriguing phase 
diagram in tilted field. 

In 1997, Zheng, Radtke, and Das Sarma (ZRD)5-^ 
predicted that the bilayer quantum Hall systems at to- 
tal filling fraction v — 2 can exhibit a novel sponta- 
neously symmetry-broken phase. ZRD performed a time- 
dependent Hartree- Fock study of spin-density excitations 
of the 11 — 2 bilayers and found that, under experimen- 
tally attainable conditions (at finite Zeeman energy and 
interlayer tunneling), the spin-density mode softens, sig- 
naling a phase transition to a new, symmetry-broken, 
quantum Hall state. The novel state was found to have a 
finite magnetization, like the simple ferromagnetic state, 
(F), that occurs in the limit of large Zeeman energy; it 
also exhibited interlayer phase coherence, like the spin- 
singlet state, [S), stabilized by strong interlayer tunnel- 
ing. In addition, the new state was found to possess anti- 
ferromagnetic correlations. ZRD dubbed the novel state 
"canted" , since schematically the state can be viewed 
as a one in which the spins in the opposite layers are 
canted away from the magnetic field in opposite direc- 
tions. The canted ground state breaks the C/(l) symme- 
try associated with rotations around the direction of the 
magnetic field; this spontaneously broken symmetry is 
behind the formation of the soft (Goldstone) spin-density 
mode found by ZRD. 

The canted state is a pure many-body state, stabilized 
by Coulomb interactions. In the absence of the inter- 
actions, there would be a first-order phase transition be- 
tween the ferromagnetic phase and the spin-singlet state. 
Coulomb interactions effectively mix the ferromagnetic 
state and the spin-singlet state around the would-be first- 
order phase transition (when the energy splitting between 
these states is small), giving rise to the canted phase. 
Because Coulomb interactions somewhat favor the ferro- 
magnetic state, a finite amount of tunneling is crucial for 
the stability of the canted phase. 

Brey, Demler, and Das Sarma (BDD)S, however, 
showed, that the canted phase can be extended to the 
region of finite Zeeman energy and infinitely small tun- 
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neling by creating a charge imbalance between the lay- 
ers. The charge imbalance can be induced by an ex- 
ternal bias voltage applied perpendicularly to the sys- 
tem. In the charge-unbalanced regime, another many- 
body phase with a spontaneously broken U(l) symme- 
try, but with vanishing antifcrromagnetic order parame- 
ter, was discovered at zero tunneling. This phase — the 
/-phase — is continuously (without a phase transition) 
connected to the antiferromangetic canted phase. Mac- 
Donald, Rajaraman, and Jungwirth (MRJ)^" pointed 
out in their thorough Hartree-Fock study of the v — 1 
bilayer phase diagram that the /-phase is akin to the 
spontaneous interlayer phase coherent phase that oc- 
curs in the v — \ bilayers in the absence of tunneling^. 
MRJ therefore suggested that the spontaneous interlayer 
phase coherence of the /-phase may lead to interesting ef- 
fects in tilted fields, closely related to the commensurate- 
incommensurate transition of the v — \ bilayers^iii^. 

In a recent paper, Burkov and MacDonald (BM)ii ex- 
plore this possibility. Indeed they find that a tilted field 
applied to a charge-unbalanced j/ = 2 bilayer system in- 
duces a new quantum phase transition within the canted 
phase. Their "unrestricted Hartree-Fock" analysis shows 
that the phase transition is first-order and is between two 
commensurate phases: At low tilt angles the phase is a 
simple canted commensurate phase, in which the layer 
degree of freedom (isospin) is commensurate with the in- 
plane component of the magnetic field. At higher tilt 
angles a phase transition occurs to a phase, in which not 
only the isospin, but also the syin becomes commensurate 
with the in-plane field (Fig.^. Thus, BM conclude, the 
phase transition is not a commensurate-incommensurate 
transition. 

In this paper, their results are extended and given a 
physical explanation. Using the "unrestricted Hartree- 
Fock" approximation, we find that the new first-order 
phase transition terminates with a quantum critical tran- 
sition embedded within the canted phase (Fig. OJ. The 
simple commensurate and the spin-isospin commensu- 
rate canted phases are continuously connected to each 
other — much like the familiar example of water and its 
vapor, the two phases possess the same symmetry prop- 
erties (the spontaneously broken U{1) in our case). To 
illustrate the novel first-order transition and the quantum 
critical transition that terminates it, we present a series of 
phase diagrams (Fig.|2J). Each phase diagram is obtained 
for fixed perpendicular-field Zeeman energy and tunnel- 
ing strength, and the axes on the phase diagrams are the 
in-plane field and the bias voltage. We find this choice 
of axes particularly suitable, since current experimental 
techniques allow to vary the bias voltage and the in-plane 
field in situ over a wide range of values. All the phases 
on a given phase diagram can therefore be accessed on a 
single sample; this should facilitate the detection of the 
novel phase transitions and the new phases. In an up- 
coming publication^^, we obtain the collective modes in 
various parts of the phase diagram, thus providing signa- 
tures of different phases and phase transitions for possible 



light-scattering experiments. 

In addition to our new Hartree-Fock results, we present 
a detailed physical discussion of the surprising behavior 
of the 1^ — 2 bilayers in tilted magnetic field. We show 
that the exotic spin-isospin commensurate canted phase 
is closely related to the /-phase, as illustrated by a com- 
parison of order parameters in the two phases (Fig. [SJ. 
In fact, the spin-isospin commensurate phase rapidly con- 
verges to the /-phase as the tilt angle is increased. We 
use the symmetry properties of the /-state to give a sim- 
ple explanation to the spin-commensuration at high in- 
plane fields. The similarities between the novel first-order 
transition of the v = 2 bilayers and the commensurate- 
incommensurate transition of the v = 1 bilayers help 
us understand many aspects of our and BM's numeri- 
cal findings, such as the confinement of the novel phase 
transition to the canted phase and the absence of the new 
transition in charge-balanced v = 2 bilayers. 

This paper is organized as follows: 

In Section m the bilayer Hamiltonian in under the 
infiuence of tilted magnetic fields is presented. 

In Section IIIII the Hartree-Fock procedure, used by 
Das Sarma, Sachdev, and Zheng (DSZ)^ to obtain the 
phase diagram of the charge-balanced v = 2 bilayer sys- 
tem in perpendicular field, is extended to the present case 
of charge-unbalanced 1^ — 2 bilayer in tilted field. A sim- 
ilar procedure was outlined by BM^^, who dubbed it the 
"unrestricted" Hartree-Fock approximation. In Sec. IIIII 
a detailed presentation of this, unrestricted, Hartree-Fock 
procedure is given. 

In Section Hvl the global phase diagram of the ly — 2 
bilayers in tilted field is presented. The phase diagram 
of the V = 2 bilayers in perpendicular field has been 
previously reported and discussed by several groups of 
authors2*2i2ii£. MRjiS in particular, have obtained the 
full Hartree-Fock phase diagrams of the = 2 bilay- 
ers for various combinations of parameters. MRJ used 
a very elucidating reduced Hartree-Fock solution, which 
elegantly captures the physics of the i/ = 2 bilayers. The 
only shortcoming of their approximation is that it does 
not always give the exact Hartree-Fock ground state that 
is crucial for obtaining the correct collective mode dis- 
persions from time-dependent Hartree-Fock (as will be 
presented in Ref . . In particular, it gives an approxi- 
mate ground state for the canted phase C, in which, as 
will be shown in Sec. [3 the most interesting phenomena 
happen when the magnetic field is tilted. In addition, 
some aspects of the phase diagram, such as the proper- 
ties of the zero-tunneling /-phase have been essentially 
uninvestigated. In Sec. IIVI some of these properties are 
discussed in detail. 

In Section the main results of this paper are pre- 
sented. In the first part of this section, the physics of the 
commensurate-incommensurate transition of the v = 1 
bilayers in tilted fields is reviewed. In the next part, 
the possibility of a similar transition in v ^ 2 is dis- 
cussed. Next, the quantitative Hartree-Fock results are 
presented and explained. The results are summarized in 



3 



Figs. 121 EI and ]E\ Figure ^ illustrates the commensura- 
tion of spin with the in-plane field at large tilt angles. 
Figure El presents a set of global phase diagrams of the 
1/ = 2 bilayers in tilted magnetic fields, which show the 
emergence and the evolution of the novel phase transi- 
tion, induced by the in-plane field. Figure [3 illustrates 
the close relationship between the novel commensurate- 
commensurate transition (with the involvement of the 
spin) and a "naive" commensurate-incommensurate tran- 
sition (in which the spin is not involved). 

Sec. IVII gives a short summary of the results presented 
in this paper. 



II. THE BILAYER HAMILTONIAN 

We model the v = 2 bilayer quantum Hall system in 
tilted magnetic field by a simple Hamiltonian which in- 
cludes the five most important aspects of the system: 
Landau level quantization, Zeeman energy, tunneling 
between the layers, an external bias voltage, and the 
Coulomb interactions between electrons. Disorder will 
be completely neglected throughout this work. 

We choose to work in the gauge M^) — 
(0, B±x, — i?||a;), where B± is the component of the mag- 
netic field perpendicular to the plane of the sample, and 
is the in-plane component of the field (the total field 

is Btotai — ^J^i_ + ^f\)- If the layers are assumed to be 

infinitely thin (an approximation used throughout this 
work), the lowest-Landau-level single-electron wavefunc- 
tions for each layer in the gauge A(r) are 
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where X ~ —kyP are the guiding centers of these 
Landau-gauge wavefunctions; I = y/hc/eB± is the mag- 
netic length, and Ly is the length of the system in the 
y-direction. Throughout this paper, we assume that the 
cyclotron energy is much larger than all other energy 
scales in the system and restrict our arguments to the 
lowest Landau level. 

In addition to the orbital degrees of freedom, X, elec- 
trons in the bilayer systems possess a spin and can be 
localized in either layer. The layer index serves as an 
additional discrete degree of freedom — the isospin — in 
bilayer systems. As a result, the electron creation opera- 
tors cj^sx labeled by three indices: the orbital index 
X, the spin index s —1 (or -|-1) and J, (or —1), and the 
layer index fi that can take on values R (or -1-1) and L (or 
— 1). To facilitate the analogy between the spin and the 
isospin, we define the electron spin and isospin operators, 
Sx and Ix- 



fj. ss' 
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where a and r are sets of Pauli matrices, and the spin 
operator is defined here in the usual mannei^. The total 
spin and isospin of the system are defined by the opera- 
tors = Y.xOx, where = I,S. 

The term in the Hamiltonian representing the Zeeman 
energy is simply 

1 
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Z'JX 



(4) 



where = g^J'B Btotai is the Zeeman splitting, and Ni 
and iV| are the total numbers of down and up spins in 
the system. Very similarly, the term representing the 
bias voltage (i.e. the difference in electrostatic potential 
between the layers) is written as 
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v^x 



X 



Nl), (5) 



where Ay is the potential difference between the lay- 
ers, and Nl and Nn are the total number of particles 
in the left and right wells respectively. The bias voltage 
clearly acts as an effective external field that couples to 
the isospin. 

Within the tight-binding approximation^"*, tunneling 
between the layers can also be expressed as an effective 
external field coupling to the isospin. When the magnetic 
field applied to the system is perpendicular to it (so that 
the electrons tunnel along the direction of the magnetic 
field and never see any magnetic flux), the tunneling term 
can be written as 
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(6) 



(cLxCLsX + cLxCflsx) , (7) 



i.e. it acts as an external field acting on the isospin in 
the /^-direction. The coefficient A^^^ is the symmetric- 
antisymmetric gap induced by the tunneling in the ab- 
sence of the other external fields and interactions. 

When the magnetic field is tilted with respect to the 
normal to the sample, the electrons tunneling between 
the layers pick up an Aharonov-Bohm phase [| / A • dl]. 
As a result, the tunneling term acquires phase factors: 

Ht = -A5AsE[e''^"''4.xCL.x + 
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^SAS 



+ e-^«"^cL^c«.x] (8) 
^[e«Qii^/++e-^'3ii^/-] (9) 
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where the operators — ± il^ are the isospin raising 
and lowering operators. The wavevector Qn is defined as 
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with d the distance between the two layersi^ (and the 
in-plane magnetic field, in our gauge, pointing in the y- 
direction). Note that |Q||| is proportional to t&nO = 
B^/B±, where 6 is the tilt angle. Interference between 
the electrons tunneling in the presence of an in-plane 
field results in the reduction of the effective symmetric 
antisymmetric gapi^ 

AsAS^A^Ase-"^^^''^'- (11) 

We note that, had we chosen a different gauge, the Gaus- 
sian decay of As^g would have remained although the 
form of this term would have been different, and the 
phase factors could disappear from the tunneling term 
and reappear in other terms in the Hamiltonian (see 

sec.rra. 

The three terms of the Hamiltonian {Hz, Hy and Ht) 
given above all describe single electrons and comprise the 
non-interacting part of the Hamiltonian 

iJo = Hz + Hv + Ht 

X 

+ ^|i5(g.Q|,X^+^g-.Q,|X^^)]^ (12) 

The Coulomb interactions between the electrons are 
taken into account by an additional term 

where intra- and mferlayer Coulomb interactions are 

= Vnd<l)^'^-^e-^\ (14) 

eq eq 

respectively, d is the distance between the layers, and 
is the area of the sample. The total Hamiltonian is 
therefore simply 

H = Ho + Hi. (15) 

III. THE UNRESTRICTED HARTREE-FOCK 
APPROXIMATION 

A. Trial ground state 

The Coulomb-interacting Hamiltonian in Eq. 1)15(1 is 
not tractable exactly, and we solve it using the Hartree- 
Fock approximation. In the usual manner, we assume 
that the many-body ground state, \G), is a Slater de- 
terminant of single-particle states and perform a func- 
tional minimization of the expectation value {G\H\G) 



with respect to these single-particle states. Under the 
assumption of translational invariance in the y-direction, 
the most general single- particle state is a superposition 
of all the combinations of spin and isospin degrees of free- 
dom — R 1 , R I, L 1 , and L I — which can be described 
as a normalized 4-spinor 

W = {wB.'f,WRi,WLT,'WLi), (16) 

Such a normalized complex four dimensional vector 
transforms under U{A). To make a v = 2 Slater determi- 
nant state, we occupy each orbital X by two particlesii 

|G)=n/ix/2VlO), (17) 

X 

where 

flx=J2<'rxcUx (18) 

creates a particle described by Eq. Hlt)|) . The requirement 
that the operators obey the fcrmionic anticommuta- 
tion relations 

{fmX, fix'} = SmnSxX' (19) 

is equivalent to an orthonormality constraint on the 

Y.iK'^xTK.x = 6n^- (20) 

Each element of C/(4) specifies four orthonormal W" 
spinors. When the filling fraction is v — 2, the two states 
with lowest mean-field energies (which we label 1 and 2 
throughout the paper) are occupied^^. 

The unrestricted Hartree-Fock ground state given by 
Eq. H17|l is very general. The freedom provided by the 
coefficients w^^^ includes a possibility of non-uniform 
states, such as stripe statesi^. However, in this paper, we 
restrict our attention to states with uniform density and 
uniform spin density, ignoring the possibility of charge 
density waves or spin density waves. Our choice is sup- 
ported by the absence of instabilities in the excitation 
spectra above the uniform-density states, which are cal- 
culated in Ref. 12. 

In order to have a state of uniform spin density and 
uniform real density in each layer, the occupation number 
Nfj,ax of each state ^aX 

N^aX = \wl,x\'+Kax\' (21) 

has to be independent of the position X . Thus, we may 
write24 

K-^x = e^^^-^^'-'z-^, (22) 

where z^^ 's are independent of position and the only po- 
sitional dependence arises in the phases. In the case of 
zero in-plane field (i3|| = Q\\ = 0), it will be very easy to 
see below that the lowest energy solution should have no 
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positional dependence of the phases, so that (p^^^X) can 
be taken to be zero. However, in the more general case of 
nonzero in-plane field a nontrivial positional dependence 
will be favored. 

Following BMii, we use a simple trial form for the 
positional phase dependence 

(j)^.{X) = Q^^X (23) 

or, equivalently, 

/nX=E(^;<^)*e"''^^"''^M<.^- (24) 

A ground state (Eq. ((171) ') with nonzero Q^^r possesses 
spin-isospin-wave order, discussed at length in Sec. IV Dl 
As was mentioned above, the proposed ground state 
(Eq. H24I) and (|17|) ) is not the most general Slater de- 
terminant (Hartree-Fock) state. Indeed, having made a 
specific choice (Eq. (j2Sl) of the positional dependence 
of the phase, this trial state is not even the most gen- 
eral state with uniform density and spin density in each 
layer. However, our analysis of the collective modes^^ 
around the ground states obtained by the minimization 
of {G\H\G) indicate the stability of these states against 
second-order transitions that cannot be described within 
the Hilbert space defined by our ansatz. The possibil- 
ity of phase transitions into a soliton-lattice state will be 
discussed in the concluding section FVll 

B. Hartree-Fock minimization 

We minimize the expectation value of the Hamiltonian 
(i.e. the trial ground state energy) with respect to the 
variational parameters z^^., and Q^a. Since 

n = l,2 

the expectation value of the ground state energy per unit 
flux in our non-uniform ansatz is 

\{G\H\G) = (26) 

fiuss' 

+ AsAsSss'[cos{{Q\i - Qrs + Qls)X)t^^ + 

+ sin((Q|| - Qn, + Qls)X) r^«,]} ^ (z^J^z^f,, + 

n=l,2 

fj.s ji=l,2 s',m=l,2 

E -^/^"^(-(Qps - Qi^sOfc) X 

fl.L^,Ss' 

m,n— 1,2 



where g is the Landau level degeneracy, g = r^/27r/^, and 
= / §e-'=^'^/V,.(fc)fcJo(M') (27) 

HM) = i(^-(^)-^-(^)) = ^^J^^8) 

so that = -ff-(O) = 1^^. The functions are 
monotonically decreasing functions of q and therefore the 
higher the wavevector of the spin-isospin-wave order, the 
lower the contribution of exchange to the ground-state 
energy. The exchange term thus favors a uniform state. 
The tunneling term, however, does not contribute to the 
ground state energy, unless Qrs — Qls = Q\\, and there- 
fore favors isospin-wave order. It is this competition that 
gives rise to the novel first-order transition discovered by 
BM^l and is presented in more detail in the next sec- 
tion. Using the results of BM"-\ we make a simplifying 
assumption that — Qlo- ^ Qi , for both ct =t and |, 
and that Q^j — — Qs, for ^ — R and L. In other 
words, we reduce the number of phase parameters from 
three to two and write 

= |g/ + (29) 

where a finite Qi indicates the presence of an isospin- 
wave order, while a finite Qs reflects the real spin-wave 
order. 

To find a variational minimum of this Hamiltonian, we 
differentiate ^{G\H\G) with respect to (zJJ^) — under 
the orthonormality constraints on (Eq. H2U|) '). The 
resulting set of minimization conditions can be arranged 
in the form of a Schrodinger equation: 

MZ" = e„Z" (30) 

where = (^aT' ^Rli ^l-], -^Li), and A/ is a 4 x 4 matrix, 
which is just the mean-field single-particle Hartree-Fock 
Hamiltonian: 

Mys>.f_,s = —Az6f^i,a^^, — AySss'T^j^ — 

- AsAsSss' [- cos((Q|| - Qi)X) r;, + 

^ X 

-f i^sin((Q||-Q,)X)T^^J + 
^ X 

fis s',m— 1,2 

X E i^'^s'T^'Hs- (31) 

n=l,2 

The Schrodinger equation (Eq. (|3()|) ) is solved 
iteratively^. At each iteration, the two eigenstates corre- 
sponding to the lowest eigenvalues are filled (i.e. chosen 



to be the states 1 and 2) . These lowest-energy eigenstates 
and are then used to obtain the matrix M for the 
next iteration. The procedure is repeated until a self- 
consistent solution is achieved. This solution — a set of 
eigenspinors Z" sorted according to their eigenvalues — 
defines the lowest energy trial state among the Slater de- 
terminants defined by Eqs. H17|l and (|24|l subject to fixed 
values of the Q^s 's. The eigenvalues e„ give the binding 
energy of a particle in the subband n, i.e. the energy lost 
when the particle is taken out of the system. The sum 
of individual binding energies does not give the ground- 
state energy; the ground state energy is evaluated using 
Eq. H27I) . The minimization of the energy of the ground 
state over the Q^s's is done last^^. Thus, we find the 
Hartree-Fock ground state in two steps: First, for fixed 
values of Qi and Qg, we minimize the expectation value 
of the ground state energy with respect to (z^^^Y . Then, 
we minimize the ground state energy with respect to Qi 
and Qg. The results are summarized in the following 
sections. 



IV. GLOBAL PHASE DIAGRAM OF THE v = 2 
BILAYERS IN ZERO IN-PLANE FIELD 

We start by considering the zero in-plane field case. 
As was mentioned in the introductory section, the phase 
diagram of the v = 1 bilayers in perpendicular field has 
been studied in some detaitf'S'^iiii. In this section, we 
highlight the properties of the zero in-plane field phase 
diagram that help elucidate the physics of the v — 2 
bilayers in tilted field and discuss a number of issues, to 
which little attention has been paid to date. 

The Zeeman splitting for a typical bilayer GaAs 
sample^ with filling fraction = 2 in perpendicular field 
is ~ O.OI(e^/eZ); it is set by the properties of the 
material (the g-factor and the effective electron mass) 
and the density of the 2DEGs. We assume that, in the 
absence of an external bias voltage, the density of the 
2DEGs in the two layers is the same. Applying a fi- 
nite bias voltage. Ay, perpendicularly to the layers, one 
can induce a charge imbalance between the layers. We 
assume that the charge imbalance can be created while 
keeping both the filling fraction and the magnetic field 
constant. The tunneling strength, A^^g is assumed to 
be unaffected by the induced charge imbalance. 

It is apparent from Eq. (|27|l that, when the magnetic 
field is oriented perpendicularly to the plane of the bilayer 
sample, i.e. =0, spin- or isospin-wave order is not 
favored. Throughout this section, we therefore assume 
that Qi = Qs = 0. 

Under these conditions, we obtain the global phase di- 
agram for the V = 2 bilayers. A cross-section of the phase 
diagram for A^ = O.OI(e^/eZ) is presented in Fig.@.. The 
phase diagram exhibits four phases: ferromagnetic 
spin-singlet S*, canted C, and the spin-isospin-entangled 
non-canted phase / (represented by the thick line along 
Ag^g = 0). As expected, when both the bias voltage 
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FIG. 1: Global phase diagram for a. v = 2 bilayer sample in 
perpendicular magnetic field. The Zeeman energy is ~ 
0.01(e'^/eZ), the interlayer spacing m d = I. The tunneling 
strength A5^3 and the bias voltage Ay are given in units 
of (e^/el). The phase S is the spin-singlet phase, F - the 
ferromagnetic phase, and C - the canted phase. The thick 
line along A^^g = represents the I phase. 



Ay and the tunneling strength A^^^ are small, the Zee- 
man energy dominates and gives rise to the ferromag- 
netic phase. In the opposite limit, of large A^^^ and 
Ay, the spin-singlet phase is stabilized. In the interme- 
diate regime, the Coulomb interactions give rise to the 
canted phase when A^^^ 7^ 0, and to the /-phase when 

The topology of the phase diagram is the same for 
all other finite values of A^. For larger values of A^, 
the phase space volume of the ferromagnetic phase in- 
creases, and the spin-singlet phase is shifted to higher 
values of A^^^ and Ay; the width of the canted phase 
(slowly) decreases. For smaller values of A^, the oppo- 
site effect takes place: the volume of the ferromagnetic 
phase decreases and the canted phase becomes wider. 
The ferromagnetic phase does not disappear from the 
phase diagram until A^ — 0, when, within the Hartree- 
Fock approximation, a many-body phase (the A^ — > 
limit of the canted phase) fills the low-A^^^-Ay region 
completelyi^^ In real bilayer samples A^ is always finite, 
while < Ag^5 < 0.08(e^/e/). As was proposed by 
Brey, Demler, and Das Sarma2-, and is clear form Fig. ^ 
by sweeping the external bias voltage one can probe the 
three phases of the u = 2 bilayer system that occur in 
the presence of finite tunneling - ferromagnetic, canted, 
and spin-singlet. In tilted fields, the /-phase can also be 
attained in the limit of a large tilt angle (see Sec lVP)) . 
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A. Ferromagnetic phase 

The simplest phase in the phase diagram is the ferro- 
magnetic phase. The ferromagnetic ground state has the 
simple form \F) — Ylx ^^rix^^l]x\^) ■ effectively a 

single-particle state, which could occur in the absence of 
interactions. The i/ = 2 system in the ferromagnetic state 
can be viewed as two decoupled spin-polarized v = 1 
monolayers. The state is clearly not interlayer phase co- 
herent, and an in-plane field therefore would not affect 
it. 



B. Spin-singlet phase 

The spin-singlet state is also a single-particle state. 
It is stabilized by large A^^^ and/or Ay. The spin- 
singlet state, IS") = Ylxi^R^^mx + ^Lcl^x)i^R^Rix + 
zlc^^I^) |0), is interlayer phase coherent unless A^^g. = 



0. 



When A^^s 



0, the spin-singlet state is simply 



X '^R]X'^RIX\ 



.\0) - a, 1/ — 2 monolayer quantum 
Hall state. In the presence of interlayer tunneling, the 
v = 2 bilayer system in the spin-singlet state can be 
viewed as two oppositely spin-polarized v = 1 bilayer sys- 
tems that possess interlayer phase coherence. The main 
difference between the = 2 bilayer system in the spin- 
singlet state and a set of two 1^=1 bilayer systems is 
that in the latter case the interactions play an impor- 
tant role alongside tunneling in creating the interlayer 
phase coherence; in the spin-singlet state of the f = 2 
bilayers interactions are not important. The phase space 
region where the interactions play an active role in = 2 
bilayers is the region of stability of the canted phase. 
The boundaries of the canted phase therefore mark the 
boundaries of the influence of the interactions. Moreover, 
unlike in the v — 1 bilayers where tunneling and exchange 
work cooperatively, in the v = 2 bilayers tunneling and 
exchange are in competition, since the Coulomb interac- 
tions favor the ferromagnetic state (due to their intra- 
interlayei anisotropy). 



C. Canted phase 

While the ferromagnetic and the spin-singlet phases 
are essentially single-particle phases, stabilized by single- 
particle fields, the canted phase is a many-body phase 
stabilized by the interactions. As was mentioned in the 
introduction, in the absence of the interactions there 
would be a first-order phase transition between the spin- 
singlet and the ferromagnetic phases. The interactions 
can lower the energy of the system by creating the 
canted phase (Fig. and the ground state energy and 
width of the canted phase depend on the strength of the 
interactions^^. As it will be shown in Sec. the impor- 
tance of the interactions in the canted phase implies that 



the phase can be non-trivially affected by the in-plane 
component of the magnetic field. 

In the canted phase, the interactions effectively mix the 
ferromagnetic and the spin-singlet states-^ giving rise to 
a finite magnetization (S^) and antiferromagnetic spin 
correlations 1(51, x Sr) — {SI x Sl)\ 0, where is 
the unit vector in the spin-up direction in the layer 
(One often defines an antiferromagnetic order parame- 
ter O^J^, where Oo^p = (5" (g) J'^), which is finite 
only in the canted phase). The U{1) symmetry associ- 
ated with rotations around is spontaneously broken. 
The canted ground state, |C), has the most general form, 
Eq. H17|l . and is spin-isospin entangled — i.e. it cannot 
be decoupled into two independent spin or isospin chan- 
nels like the ferromagnetic or the spin-singlet ground 
states. The canted ground state is therefore interlayer 
phase coherent, which is confirmed by a non-zero value 

of ((/-)2 + (/!')2). 



D. /-phase 

While the canted phase has attracted the most 
attention£i2i2ii2ii&, the many-body phase that occurs in 
charge-unbalanced systems in the absence of interlayer 
tunneling is no less interesting. In the absence of tun- 
neling, the interactions give rise to a spontaneously inter- 
layer phase coherent /-state^-^" (with the antiferromag- 
netic order parameter, Ozx = 0). The interlayer phase 
coherence of the /-state is spontaneous since the single- 
particle fields A^ and Ay do not stabilize interlayer 
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FIG. 2: Energy profile of charge-balanced u — 2 bilayers with 
Zeeman energy Az ~ 0.05(e^/ei), interlayer spacing d — I, 
and a range of tunneling amplitudes. The dashed line repre- 
sents the energy of the ferromagnetic state. The long dashed 
line is the energy of the spin-singlet state. In the absence 
of interactions, there would be a first order phase transition 
at the intersection of the two dashed lines. The interactions 
effectively "smooth out" the profile by stabilizing the canted 
phase. 
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phase coherent states. The /-ground state has a simple 

form |/) = rix ^^RT xi^Rl'^^Rlx ~^ '^L'\c'L-tx)\^) ' which man 
ifests the interlayer phase coherence or the state, its spin- 
isospin entanglement, and the U{1) spontaneous symme- 
try breaking. Simply speaking, the U{1) symmetry of 
the /-state is the freedom to choose the relative phase 
of zj^i and z^j in the expression for |/). One can also 
formally consider the U{1) symmetry associated with ro- 
tation around /^ or around (applying e*^^' or e*^'^" 
clearly gives the desired effect). It is important to note 
that the state |/) doe not break the U{1) x U{1) sym- 
metry of spin and isospin rotations completelySiiS. This 
state is an eigenstate of /^ -I- with an eigenvalue g 
(the degeneracy of the Landau Level), but it mixes states 
with different /^ — quantum numbers. So the state |/) 
breaks the U{1) x U(l) symmetry down to the diagonal 
[/(I). 

The feature of the state |/) to be an eigenstate of /^ -I- 
will prove useful in the following discussion of the 
many-body phases in the presense of an in-plane field. In 
fact, all the possible ground states of the 1^ — 2 bilayers in 
the absence of tunneling are eigenstates of + S^. If the 
bias voltage and the Zeeman energy are positive, all the 
possible zero-tunneling ground states are (locally — for 
a given orbital quantum number X) linear combinations 
of a = +1 spin triplet and an /^ = -1-1 isospin triplet 
(the ferromagnetic and the spin-singlet states being the 
two extremes with only the spin triplet or isospin triplet 
contributing, respectively). In this case, all the zero- 
tunneling ground states are eigenstates of /^ -|- 5*^ with 
an eigenvalue g, where g is the Landau level degeneracy. 

Even though I^+S^ is a formal construction, whose ex- 
pectation value cannot be directly measured, it can help 
advance our understanding of the physics of the /-phase 
and of the canted phase in the presence of small interlayer 
tunneling. Thus, since applying a rotation e**'^ ^ to 
the state |/) generates a trivial phase factor, /^ + can 
be treated as the effective direction in which the /-state 
"points" . In contrast, if tunneling is added to the system, 
I^ + S^ does not commute with the Hamiltonian. The op- 
erator e*^*^^ "'"■^ ^ then generates a non-trivial rotation of 
the canted state around P + S^. When interlayer tunnel- 
ing is small, the canted ground state has a large overlap 
with an /-state, and can be visualized in its Hilbert space 
as pointing slightly away from the /^ + 5^ direction. The 
relevance of this picture to the v — 2 bilayer physics will 
become more clear in the next section fSec lVD)) . 

The relation between the canted phase and the /-phase 
is reminiscent of the relation between phases of the i> — 1 
bilayers in the presence and absence of tunneling. The 
main similarity is that, when there is no tunneling be- 
tween the layers, both the i' — 1 and the v = 2 bilayers 
support a spontaneously interlayer phase coherent phase. 
There are, however, marked differences. The main differ- 
ence is that the symmetry properties of the canted phase 
and the /-phase are the same, while in the v = \ bilayers 
spontaneous symmetry-breaking occurs only in the ab- 
sence of tunneling. Moreover, unlike the zero-tunneling 



phase of the v = \ bilayers, the /-phase is spin-isospin 
entangled. As we will show in the next section, be- 
cause of their similar nature, both the v — \ and v — 2 
systems undergo essentially similar phase transitions in 
tilted magnetic field; the systems' differences, however, 
lead to surprisingly different behavior of the v — \ and 
V — 2 bilayers in the vicinity of the phase transitions. 



V. GLOBAL PHASE DIAGRAM OF THE v = 2 
BILAYERS IN TILTED FIELD 

As was discussed in the previous section, the physics of 
the V = 2 bilayers is very rich. The v = 2 bilayers exhibit 
a host of many-body phenomena, such as spontaneous 
symmetry breaking, spontaneous interlayer phase coher- 
ence, and spin-isospin entanglement. To further explore 
the many-body nature of the phases and phase transi- 
tions of the V = 2 bilayer system, in this section we 
study the behavior of the system in the presence of a 
tilted magnetic field. 

Tilted magnetic fields have been successfully used to 
study the role of both spin and layer degrees of freedom 
in quantum Hall physics. In thin monolayer systems, the 
tilted field technique was used to investigate the spin- 
unpolarized fractional quantum Hall ground states. The 
technique is based on the fact that, in the infinitely thin 
limit, the orbital motion in a 2DEG depends solely on 
the perpendicular component of the magnetic field, /3_l, 
while the Zeeman energy is proportional to the total field, 
B. The Zeeman energy can therefore be increased in- 
dependently of the effective interactions in the quantum 
Hall monolayer by adding an in-plane field B\\. In bilayer 
systems, the presence of an in-plane field affects not only 
the Zeeman energy but also the relative orbital motion 
in the two layers and, therefore, the tunneling between 
the layers and interlayer interactions. 



A. Commensurate-incommensurate transition in 
the V = \ bilayers: An overview 

The coupling of the tilted field to the layer degree of 
freedom is easier to demonstrate using the thoroughly 
studied v — \ bilayersSiSi^i^ as an example. In i/ = 1 bi- 
layers the spin degree of freedom is frozen out by the fer- 
romagnetic exchange^ and, in the absence of an external 
bias voltage, the physics of the system is fully determined 
by the interplay of tunneling and Coulomb interactions. 
As was discussed above, tunneling (Eq. lO) can be con- 
sidered as an external field that couples to /^; Coulomb 
interactions give rise to a charging energy (from the direct 
term) and an anisotropic isospin-stiffness, both of which 
favor the isospin- a;y planed. When the layers are sepa- 
rated by a distance comparable to the distance between 
the electrons in a single layer, the anisotropic Coulomb 
interactions support a spontaneously interlayer phase co- 
herent ground state even in the absence of tunneling^. 
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The U{1) symmetry associated with rotations around P 
is spontaneously broken. TunneHng, always present to 
some degree in real samples, breaks this symmetry, but 
it does not destroy the interlayer phase coherence. In- 
stead, in the absence of an in-plane field, tunneling acts 
cooperatively with the interactions to stabilize the inter- 
layer phase coherent state^ii.. 

A finite in-plane field introduces a competition be- 
tween tunneling and the Coulomb interactions in ly = 1 
bilayers^. In the presence of an in-plane field, the tun- 
neling term (Eq. lO) now favors an isospin-wave ground 
state, in which the isospin twists around the /^-direction 
with the wavevector (5||. Exchange interactions, on the 
other hand, favor a uniform configuration. The com- 
petition between tunneling and exchange results in the 
commensurate-incommensurate transition between the 
isospin-wave commensurate and the uniform incommen- 
surate phases. 

The picture of the commensurate-incommensurate 
transition outlined above is somewhat simplistio^. The 
transition, in fact, does not happen directly from the 
commensurate to the incommensurate state, but in- 
stead occurs through formation of a soliton-lattice phase. 
When, as the in-plane field is increased, the losses 
in exchange energy due to twisting become approxi- 
mately equal to the expectation value of the tunnel- 
ing term, the system can optimize its ground state en- 
ergy by forming a soliton in an otherwise commensurate 
sta.te^'^'i^^. By forming a soliton the system recovers 
some exchange energy while saving most of the tunneling 
energy. The commensurate-incommensurate transition 
by means of soliton-formation is the Talapov-Pokrovsky 
commensurate-incommensurate transition^ . 

The microscopic description of the physics of a bi- 
layer system in an in-plane field depends on a particu- 
lar choice of the gauge, but the underlying picture of an 
induced competition between tunneling and exchange is 
valid for any gaug^. Thus, for example, in the gauge 
A = {0,B±x — i3||Z,0), the tunneling electrons acquire 
no Aharonov-Bohm phase and the tunneling term in the 
microscopic Hamiltonian is unaffected'^' However, in 
this gauge it is the interlayer exchange that acquires an 
additional phase in this case, so that it is the interaction 
term that now favors an isospin-wave. 



B. Possibility of commensurate-incommensurate 
transition in = 2 bilayers 

To create a framework for the interpretation of our 
numerical results, we start by a qualitative discussion of 
the behavior of the 1^ — 2 bilayers in tilted field. We 
consider the theoretical possibility of the commensurate- 
incommensurate transition in the i> — 2 bilayers, much 
like it has been done for the v — 1 bilayers: We assume 
for simplicity that the commensurate-incommensurate 
transition happens between a commensurate state that 
maximizes tunneling at the expense of exchange (i.e. 



Qi = Q\\ and Qs = 0) and an incommensurate state 
that maximizes exchange while losing all the tunnel- 
ing energy (Q/ = Qs — 0). This is the scenario 
that we called simplistic in our discussion of the v — 
1 bilayers since better mathematical models indicate 
that the commensurate-incommensurate transition oc- 
curs through a succession of soliton lattice phases rather 
than directly. Nevertheless, the "naive" commensurate- 
incommensurate transition serves as a good simple ap- 
proximation that allows one to make experimentally rel- 
evant predictions and explanations without any numer- 
ical calculations. We therefore expect that the naive 
commensurate-incommensurate approximation will cap- 
ture some of the physics of the ly = 2 bilayers in tilted 
fields as well. 

The naive commensurate-incommensurate transition 
can happen only in the canted phase, since it is the only 
phase in which tunneling and exchange are comparable 
to each other. The other two phases — the ferromag- 
netic phase and spin-singlet phase — essentially mark 
the regions in the phase space where interactions are less 
important than the single particle fields. We therefore ex- 
pect the spin-singlet phase to be always commensurate. 
The ferromagnetic state is not interlayer phase coher- 
ent and therefore the in-plane field cannot affect it. The 
canted phase however can be commensurate at lower in- 
plane fields, but it may turn incommensurate when the 
in-plane fields are so high that the losses in exchange 
become greater than the contribution of the tunneling 
term. 

The naive commensurate-incommensurate transition 
however cannot happen in charge-balanced v — 2 hi- 
layers. In order for the commensurate-incommensurate 
transition to happen, there should exist an incommensu- 
rate state, which would become lower in energy than the 
corresponding commensurate state as the tilt angle is in- 
creased. Since the tunneling contribution to the energy 
of the naive incommensurate state is zero, the incommen- 
surate state is equivalent to the state of the system with 
no interlayer tunneling (all other parameters unchanged). 
In charge-balanced v — 2 bilayers, the zero-tunneling 
ground state is always ferromagnetic. As is clear from 
Fig. 12 the energy of a ferromagnetic state created in a 
system with no tunneling and given intralayer interac- 
tions is always higher than that of a canted state created 
by adding a finite amount of tunneling to this system 
(all other things being equal). Therefore, there can be 
no naive commensurate-incommensurate transition in a 
charge-balanced v — 2 bilayer system. 

The situation changes if a finite bias voltage can be 
applied to the system. As was mentioned in Sec. IIVI in 
the presence of finite bias voltage one can obtain an in- 
terlayer phase coherent phase — the /-phase — even in 
the absence of tunneling. The ground state energy of the 
/-state depends not only on the relative strength of the 
Zeeman energy and the external bias voltage, but also 
on the interlayer Coulomb interactions. Therefore, it is 
straightforward to argue that, in charge- unbalanced sam- 
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pies with small tunneling amplitudes, a commensurate- 
incommensurate transition may be possible: Let us con- 
sider such a sample. In perpendicular field, the many- 
body state of the sample is a canted state. Because the 
tunneling amplitude is small, this canted state is just a 
slight perturbation of the corresponding zero-tunneling I- 
state. The interaction energies of the two states are very 
close; the energy of the canted phase is lower largely due 
to the tunneling term. As the magnetic field is tilted, the 
isospin starts twisting commensurately with the in-plane 
field, thereby losing some interlayer exchange. When the 
losses in exchange energy become equal to the tunneling 
energy, the energies of the canted commensurate state 
and the corresponding zero-tunneling /-state become ap- 
proximately equal. At slightly higher tilt angles, a transi- 
tion to the /-state clearly would occur. In anticipation of 
this discussion, at the start of the chapter, we named the 
zero-tunneling many-body phase the /-phase, in which / 
stands for "incommensurate". 

We also can expect that the in-plane field will not only 
induce a new transition, but also affect the location of the 
second-order phase transitions between the spin-singlet, 
the canted and the ferromagnetic phases in the phase 
space. The reduced effective interlayer interactions in 
the commensurate states will destabilize the canted com- 
mensurate phase, and its two phase boundaries will move 
toward each other as the in-plane field is increased. The 
relative position of the boundaries of the incommensu- 
rate phase will be nearly constant with the increasing 
tilt angle. At larger in-plane fields another effect, which 
we so far have have ignored, will become important — 
the dependence of the Zeeman energy and the tunneling 
amplitude on the in-plane field. The increasing Zeeman 
energy and the decreasing tunneling will eventually drive 
the system into a ferromagnetic state as the tilt angle be- 
comes very large. (At very large in-plane fields, however, 
our infinitely-thin layer approximation loses its validity 
completely, and therefore this discussion becomes purely 
theoretical.) 

In the next subsection we explain that the 
commensurate-incommensurate transition at = 2 is 
sufficiently different from the simple picture presented 
above. Hence, the arguments presented here should be 
considered only as a motivation for a more detailed dis- 
cussion in subsecuent sections. 



C. Global phase diagram of v = 2 bilayers in tilted 
magnetic field 

The Hartree-Fock phase diagrams are presented in 
Fig. O The axes on the phase diagrams are the bias 
voltage. Ay, and the in-plane field wavevector We 
find this choice of axes convenient, since current experi- 
mental techniques allow to tune both the bias voltage and 
the in-plane field in situ across a wide range. The other 
parameters of a bilayer sample — the perpendicular-field 
Zeeman splitting, A^, the perpendicular-field tunneling 



amplitude, A^^g, and the distance between the layers, 
d — are, to a good approximation, intrinsic to a given 
sample. For each phase diagram, we fix these parame- 
ters at values typical of real samples: A^ — 0.01(e^/eZ), 
Ag^g < 0.08(e^/eZ), and d — I. Each phase diagram 
in Fig. 13 therefore corresponds to a single sample with 
A^ = 0.01(e^/£Z), d = I, and the value of A^^g given in 
the lower right corner of the phase diagram. 

The unrestricted Hartree-Fock calculation does pro- 
vide evidence for a new phase transitionii (the dashed 
line in Fig. ISJ, which possesses the properties we 
qualitatively predicted for the naive commensurate- 
incommensurate transition: The first-order transition oc- 
curs only within the canted phase and only in the pres- 
ence of a finite bias voltage. Moreover, the canted com- 
mensurate phase shrinks as the in-plane field is increased, 
but the decrease in the width of the phase stops after 
the first-order transition. However, instead of the in- 
commensurate phase, we find an interesting, what seems 
like a doubly-commensurate phase, in which both the 
isospin and the spin components are commensurate with 
the in-plane field^-'^. That is to say, throughout the in- 
terlayer phase coherent region — in the phases SC, CI 
and C2 on Fig. 13 — the wavevector of the isospin-wave 
Qi — Q\\- In the canted phases CI and C2 however 
the spin- wave wavevector is also nonzero. It is almost 
zero in the CI phase, except near the phase transition 
boundary, but it is close to Qs = Q\\ in the C2 phase 
(Fig. 21 and Ref . ITlh . The phase transition between the 
two canted phases is first-order, terminating at a critical 
point. The onset of the first-order transition occurs at 
a critical tunneling amplitude A°^g « 0.015(e^/e/). As 
the tunneling amplitude A^^g is increased, the C1-C2 
transition becomes more prominent and a higher in-plane 
field is needed to induce it. The presence of the in-plane 
component of the magnetic field thus leads to a phase 
transition that is clearly related to the commensurate- 
incommensurate transition, but possesses some unex- 
pected properties that invite a physical explanation. 

We emphasize that a simple picture of the 
commensurate-incommensurate transtion should not be 
carried directly from = 1 to = 2. In the former case 
it appears as a result of the competition between the sin- 
gle particle tunneling energy and the exchange part of the 
Coulomb interaction. At — 2 we have Qi — Qy in all of 
the canted phase (both CI and C2), which optimizes the 
tunneling term. The origin of the C1-C2 transtions is the 
competition of the exchange terms in equation (31). The 
wavevectors of the exchange terms in this equation are 
given by Qs, Q\\ - Qs, and Q\\ + Qs (we used Qi = Q\\). 
So different exchange terms would be minimized for dif- 
ferent values of Qs- It is useful to consider the variational 
energy of the ground state in (31) as a function of Qs for 
different points in the phase diagram, E{Qs). When we 
start near the base of the first order transition and far 
from the critical point (see Fig 3), the function has two 
local minima: one for Qs close to zero and the other for 
Qs close to Qii . When we are on the CI side of the transi- 
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tion, the former is the global minimum, and when we are 
on the C2 side, the latter corresponds to the true ground 
state. As we move along the first order line toward the 
critical point (by increasing the gate voltage), the posi- 
tions of the two local minima are moving together until 
at the critical point the merge into a single minimum. 
Anywhere above the critical point the system has only 
one local minimum in E(Qs). It is also useful to point 
out that there is a region of metastability of CI and C2 
phases around the first order line separating them. We 
expect interesting hysteresis effects to occur in this re- 
gion. 



D. Canted commensurate phases in i/ — 2 bilayers 

The most surprising part of the phase diagram is the 
C2 phase, in which both the spin and the isospin de- 
grees of freedom are commensurate with the in-plane 
fielcUi. The location of the C2 phase on the phase di- 



agram, so similar to that expected of the incommensu- 
rate phase, strongly suggests that, despite its apparent 
complexity, the C2 phase is simply related to the naive 
incommensurate phase. The close relationship between 
the C2 phase and the /-phase becomes more clear if 
one considers the expectation values of spin and isospin 
operators. In Fig. we plot three expectation values: 
(/§) = Exie"'"^^ {I' + S') for 

a sample with = O.Ol(eVd), A^^^ = 0.06(e2/d), 
and the bias voltage held at Ay = 0.8{e^/el) as we move 
across the first-order transition by increasing The 
operator Iq is the tunneling operator in the Hamilto- 
nian of the system, and therefore its expectation value 
indicates if tunneling contributes to the energy of a par- 
ticular state. The expectation value (Iq) is zero in the 
naive incommensurate state. The naive incommensurate 
state, the /-state, is an eigenstate of the operator /^ -I- 
with the eigenvalue g. Thus, ^{I^ + S^) = 1 in this 
state. To distinguish the /-state from the ferromagnetic 
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FIG. 4: Evolution of the spin-wave wavevector Qs (solid 
line) as a function of the in-plane field wavevector Qn (given 
in units of l/l). The dashed line is Qi — Q\\, given for com- 
parison. The wavevector Qs is small in the phase CI (low Q||) 
and abruptly jumps to Qs ~ Q\\ at the phase transition to 
C2. The Zeeman energy in this figure is A% — 0.01(e'^/eZ), the 
interlayer spacing is d = I, the tunneling constant is A^^s — 
0.06(e^/£/), and the external bias voltage is Ay = 



O 0.5 - 












incommensurate transition. As expected, after the tran- 
sition into the incommensurate phase, ^{Iq) — ^^nd 

^{I^ + S^) — 1. The thick lines represent the expec- 
tation values that we obtain for the system allowed to 
undergo the C1-C2 transition. The expectation values 
obtained from the full unrestricted Hartree-Fock solution 
and those obtained under the assumption of the naive 
commensurate-incommensurate transition exhibit strik- 
ingly similar behavior. The main difference is that, unlike 
in the /-phase, the expectation value of the Iq operator 
in the C2 phase is small but finite. This means that 
there is a contribution from the interlayer tunneling to 
the ground state energy in the C2 phase. It is therefore 
clear that the C2 phase is the optimized version of the 
/-phase much like the soliton-lattice phase of the v = 1 
bilayers is the optimized v = 1 incommensurate state. 

Unlike the v = \ bilayer system, the v = 2 bilayer 
system possesses the spin degrees of freedom which it 
can use to optimize its ground state energy around the 
commensurate-incommensurate transition. The system 
can satisfy the Aharonov-Bohm phases in the tunneling 
term by winding either around I' or F . Winding 
around /^ -f 5^ clearly does not affect the Zeeman and 
the bias voltage terms and in some circumstances can 
cost less exchange energy: The /-state is an eigenstate 
of the /^ -I- 5^ operator and, as we argued in Sec. IIVDI 
is invariant under rotation around -\- . If the /- 
state is somehow perturbed so that it becomes slightly 
canted, for example because a small amount of tunneling 
is present, the state is no longer invariant under rota- 
tions around /^ + , but a precession around /^ -I- S^- 
does not cost much exchange energy. (A useful analogy 
is a Heisenberg ferromagnet with all spins pointing in 
the same direction, {S^) = +1/2. If the spins are made 
to tilt away from the positive S'^-direction and precess 
around it, very little spin-stiffness energy is lost.) As the 
in-plane field is increased, the winding around /^ -I- 
becomes faster and tighter. This causes the phase C2 to 
asymptotically approach /. 



FIG. 5: Expectation values (/q) (solid lines), {P) (dashed 
lines), and {P + 5*^) (dash-dotted lines), per fiux quantum, 
across the C1-C2 (thick lines) and the naive commensurate- 
incommensurate (thin lines) phase transitions. The Zeeman 
energy is A^ = Q.Q\{e^ /el), the interlayer spacing is d = Z, the 
tunneling constant is Asas = 0.06(e^/£Z), and the external 
bias voltage is Ay = 



state and the fully charge-unbalanced spin-singlet state, 
which are also eigenstates of F + 5'^, we also plot (/^), 
which satisfies < -iF) < 1 in the /-state. The thin- 
ner lines in Fig. [S] represent the expectation values ob- 
tained under the assumption that the system undergoes 
the naive commensurate-incommensurate transition (the 
spin- wave wavevector, Qs, is held at 0). The discontinu- 
ity in the expectation values marks the commensurate- 



The CI and C2 phases are both canted and have the 
same symmetry properties. It is therefore not surpris- 
ing that they are connected on the phase diagram, i.e. 
they are essentially the same phase. The qualitative dif- 
ference between CI and C2 is in the involvement of the 
spin degree of freedom in the quenching of the in-plane 
magnetic field. In our mean-field solution the wavevector 
of the spin- wave does not always jump between the quali- 
tatively understood cases of Qs = and Qs = Q\\, but it 
can change gradually. The wavevector Qs changes grad- 
ually when CI turns into C2 via a cross-over, at larger 
values of Ay, when the canted state has a large overlap 
with the /-phase (|(^^ + S^) close to 1). In fact, when 

^SAS small, ^%as — 0.015(e^/£Z), the first order 

transition disappears altogether, since the canted phase 
is so close to the /-phase that the canted phase always 
has a C2 flavor to it. 
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FIG. 6: Schmatic representation of the comensurate CI and C2 states. In the left figure, an /-state is represented. The /-state 
is shown in Sec. IIV Dl to be an eigenstate of the operator + . Thus, the / state is represented as a vector pointing in 
the + direction. In the presence of tunneling, the system can satisfy the Aharonov-Bohm phases in the tunneling term 
by winding around either the + axis, yielding the C2-commensurate state shown above, or by winding around the P 
axis yielding the Cl-commensurate state. The C2-state asymptotically approaches the /-state in large in-plane field where the 
winding becomes very fast and therefore must be very tight (due to "spin stiff'ness"). 



VI. CONCLUSIONS 

To summarize, we have obtained the global phase dia- 
gram of the = 2 bilayers in tilted magnetic field (Fig.|2I). 
We found that, in charge-unbalanced v — 2 bilayers, a 
finite in-plane component of the magnetic field can in- 
duce a first-order phase transition between two commen- 
surate canted phases CI and C2. The phase CI possesses 
isospin-wave order, commensurate with the in-plane field; 
in the phase C2 commensurate spin-wave order is in- 
duced alongside with the isospin-wave order. Both CI 
and C2 phases spontaneously break a global U{1) sym- 
metry, and are technically the same phase. Indeed, in the 
phase diagrams in Fig.O phases CI and C2 are topolog- 
ically connected, and the first-order transition between 
them terminates at a critical end-point. 

The physics of the commensurate canted phases was 
discussed in detail in this paper. The behavior of the 
1/ = 2 bilayers in tilted magnetic fields was compared to 
that of their ~ 1 counterparts. The phase CI was found 
to be analogous to the commensurate phase of v — 1 
bilayers, while the phase C2 was linked to the incom- 
mensurate phase. As was predicted by MRJiS, the U{1)- 
symmetry-broken /-phase, which had been predicted to 
exist in the absence of tunneling in charge-unbalanced 
V = 2 bilayer systems, was found to play the role of a 
"naive" incommensurate phase in v ~ 2 bilayers (akin 
to the "naive" — translationally invariant — incommen- 
surate phase of = 1 bilayers (see Sec. IV Al) ). In this 
paper, the C2 phase was argued to be an optimization 
of the naive incommensurate phase, much like the soli- 
ton lattice phase in v = 1 bilayers in tilted fields is an 
optimization of the naive incommensurate phase in this 
system. 

The rapid convergence of the spin-isospin commensu- 
rate canted phase to the /-phase can be used to study the 
/-phase. Much like the canted phase, the /-phase also 



possesses a number of intriguing many-body properties. 
However, the /-phase can occur only in the absence of 
interlayer tunneling - a condition impossible in a typical 
bilayer sample. Tilting the magnetic field allows one to 
access the /-phase in an experimental setting and study 
its properties. 

The possibility of the formation of an interim soliton 
phase around the C1-C2 first order phase transition can- 
not be ruled out with certainty in our approximation. 
However, the energy of a soliton phase in v = \ bilayers 
converges to the energy of the corresponding naive (trans- 
lationally invariant) incommensurate phaseiSiSS, much 
more rapidly than does the energy of C2 converges to the 
energy of the corresponding /-state. We may therefore 
conclude that, even if the soliton phase in u = 2 bilay- 
ers is possible, it will not occupy a significant amount of 
phase space. 

In this paper we present a heuristic argument and 
numerical evidence that no new transition occurs in 
charge-balanced v = 2 bilayers. This is indeed consis- 
tent with the inelastic light-scattering results by Pelle- 
grini et alii2i2i: In their experiments, Pellegrini et al. 
used the tilted-field technique to sweep over a range of 
Zeeman energy in situ. No perpendicular bias voltage 
was applied to the bilayer system; the maximum tilt an- 
gle was 9 = 45°. Pellegrini et al. obtained encouraging 
evidence of the existence of the expected phase transi- 
tions between the spin-singlet and the canted phases, as 
well as between canted and ferromagnetic phases. No 
other transitions have been reported. 

The Hartree-Fock approximation, which we used to 
obtain our results, had been shown to be robust 
for the u = 2 bilayers in perpendicular ficldsiS*2S. 
The phase diagrams obtained in the Hartree-Fock ap- 
proximation closely match those obtained using exact 
diagonalization^^. While the Hartree-Fock approxima- 
tion overestimates the size of the canted region on the 
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spin-singlet side, it reproduces the boundary between 
the canted and the ferromagnetic phases apparently 
exactl]^. Since the novel phase transition occurs closer 
to the ferromagnetic side of the canted phase, it is reason- 
able to assume that the quantum fluctuations, not taken 
into account in the HF approximation, will not wash it 
out. The quantum fluctuations will probably effectively 
renormalize the canted phase and make the first order- 
phase transition terminate closer to the ferromagnetic- 
canted line. 
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